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Baryon Acoustic Oscillation (BAO) surveys will be a leading method for addressing the dark 
energy challenge in the next decade. We explore in detail the effect of allowing for small amplitude 
admixtures of general isocurvature perturbations in addition to the dominant adiabatic mode. We 
find that non-adiabatic initial conditions leave the sound speed unchanged but instead excite differ- 
ent harmonics. These harmonics couple differently to Silk damping, altering the form and evolution 
of acoustic waves in the baryon-photon fluid prior to decoupling. This modifles not only the scale 
on which the sound waves imprint onto the baryon distribution, which is used as the standard ruler 
in BAO surveys, but also the shape, width and height of the BAO peak. We discuss these effects in 
detail and show how more general initial conditions impact our interpretation of cosmological data 
in dark energy studies. We find that the inclusion of these additional isocurvature modes leads to a 
decrease in the Dark Energy Task Force figure of merit (FoM) by 46% i.e., FoM/so ~ 0.54 x FoMad 
and 53% for the BOSS and Adept experiments respectively when considered in conjunction with 
Planck data. We also show that the incorrect assumption of adiabaticity has the potential to bias 
our estimates of the dark energy parameters by 2.7a {2.2a) for a single correlated isocurvature mode 
(CDM isocurvature), and up to 4.9(t {5.7a) for three correlated isocurvature modes in the case of 
the Boss (Adept) experiment. We find that the use of the large scale structure data in conjunction 
with CMB data improves our ability to measure the contributions of different modes to the initial 
conditions by as much as 95% for certain modes in the fully correlated case. 



I. INTRODUCTION 

Although the standard model of cosmology based on 
ACDM has not changed fundamentally in the last decade, 
there has been a remarkable refinement in our knowledge 
of the parameters describing the model. For example, the 
original supernova results gave only limits of flm < 1.5 
at 2(7 assuming a general ACDM model [l| while the lat- 
est results from the WiggleZ Baryon Acoustic Oscillation 
(BAO) survey, together with WMAP and Union2 super- 
nova data now give = 0.29 ± 0.04 at 2tT 

As a result of this progress it has become obvious that 
systematic errors are a key issue in pushing the frontier 
further. For example, in the case of supernovae there are 
important systematic errors related to the lightcurve fit- 
ter used which currently leads to shifts in the dark energy 
equation of state of about Aw ~ 0.1 0. There may be 
additional supernova systematics such as the existence of 
Type la subpopulations and correlations between abso- 
lute magnitude and host galaxy type (see e.g. [3|). 

BAO have their own associated systematic errors, such 
as nonlinear effects which potentially bias or shift the 
BAO peak, although these are believed to be fairly small 
and possible to calibrate through theoretical modeling 
and N-body simulations Q. However there is another 
theoretical systematic due to isocurvature perturbations 
that has recently received attention 0, 0] . Depending 
on how general one allows the primordial isocurvature 



admixture to be, there can be a significant impact on the 
ability of future BAO surveys to constrain dark energy 
even if one imposes the constraint that the isocurvature 
modes be undetectable by Planck alone Q. This shows 
that at least in the next generation of surveys one will not 
be able to decouple the search for dark energy with BAO 
from an understanding of the early universe, a subtlety 
that does not affect supernovae surveys. 

The key reason that even small correlated isocurvature 
modes cause a problem for BAO surveys is that they alter 
the way in which the BAO peak appears in the two-point 
correlation function of baryons, and hence, of galaxies. In 
the simple adiabatic model the BAO peak is controlled 
by the sound horizon, the distance that sound waves can 
propagate in the early universe from the time of inflation 
to decoupling. This characteristic scale depends only on 
the sound speed in the standard adiabatic picture: 

where Cs{z) = 1/^3 (1 + i?6/(l + zj) and Rt = 

31500 ujb {T^^-^^/2.7 K) which can be measured accu- 
rately with the Cosmic Microwave Background (CMB). 
By comparing the size of this scale at the time of decou- 
pling and its angular size at late times we can learn about 
the expansion history of the universe and measure cos- 
mic distances, and hence constrain models of dark energy 
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In [6| , we found a clear degeneracy between the impact 
of dark energy models and non-adiabatic initial condi- 
tions on the galaxy correlation function. In this paper 
we explain in depth why small amplitude but general ad- 
mixtures of correlated isocurvature modes can have such 
a strong impact on the cosmological constraints based 
on BAO surveys. Wc show that relaxing the assump- 
tion of adiabaticity and allowing fractions of isocurvature 
modes affects the development of the acoustic waves in 
the baryon-photon fluid. The isocurvature modes excite 
different harmonics which in turn, couple differently to 
Silk damping, and in so doing, modifies both the scale 
on which the sound waves imprint on the baryon distri- 
bution and the shape of the BAO peak. 

This paper is arranged as follows; in section |lll we 
study the evolution of the baryon density contrast under 
different initial conditions and how the structure of the 
BAO peak is altered. In section Hill a Fisher matrix for- 
malism is implemented in order to quantify the impact 
of these changes on the forecasted errors on the dark 
energy parameters from two BAO experiments, namely 
Boss [m and Adept As a prior, we include the 
information from the high-resolution CMB temperature 
anisotrop y a nd polarization spectra from the Planck 
Surveyor |53| , which should provide stringent constraints 
on the amount of isocurvature in the initial conditions. 
We also conduct a study of the potential bias in our esti- 
mates of the dark energy parameters that can result from 
an incorrect assumption of pure adiabatic initial condi- 
tions. Lastly, we show that constraints on the isocurva- 
ture parameters can be derived from BAO surveys. We 
discuss our conclusions in section ITVl 



II. THE BAO PEAK WITH ADIABATIC AND 
ISOCURVATURE INITIAL CONDITIONS 

The BAO peak is sensitive not only to the matter 
content of the universe, but also to the character of 
the primordial perturbations. The features of the BAO 
peak such as the location, width and amplitude are 
mainly dictated by the time evolution of the baryon 
density contrast Sb from the post-inflation period to 
photon-baryon decoupling. In turn, the time evolution 
of the baryon density contrast during the pre-decoupling 
period depends on the initial configuration of the 
primordial perturbations in the different species at the 
end of inflation. 

In the simplest scenario the perturbation affects all the 
cosmological species such that the relative ratios in the 
number densities remain unperturbed, exciting the adi- 
abatic mode (AD). Although adiabatic initial conditions 
are a natural feature of single-field inflationary models 
[3^ . it has been shown [2^ that four regular isocurva- 
ture (ISO) modes are allowed in addition to the adiabatic 
(AD) mode. These isocurvature modes are character- 
ized by variations in the particle number ratios but with 



vanishing curvature perturbation, with different isocur- 
vature modes excited depending on the species that are 
initially perturbed. These are namely the cold dark 
matter isocurvature (CI) mode, the baryon isocurvature 
(BI) mode, the neutrino isocurvature density (NID) and 
the neutrino isocurvature velocity (NIV) mode. While 
isocurvature modes are more difficult to physically mo- 
tivate, the possibility of correlated isocurvature fluctua- 
tions is allowed given current cosmological data [2^ [2^ 
and interesting to consider. 

We will show that different modes of the primordial 
perturbations excite different harmonics and these har- 
monics couple differently to the Silk damping, thereby 
altering the characteristic BAO scale at photon-baryon 
decoupling. After decoupling, baryon fluctuations on 
scales larger than the Jeans length Xj slow down in 
the rest frame of the cold dark matter (CDM), falling 
into the CDM potential wells, and eventually tracing 
the CDM, while on scales below Xj the fluctuations 
still oscillate, independently of the initial conditions. 
In order to study the features of the BAO peak for 
different modes, we consider the time evolution of the 
photon-baryon fluid in the tight-coupling regime. 

In this regime, photons and baryons are treated as per- 
fect fluids. The subscripts 5, c, 7 and v respectively de- 
note the baryons, CDM, photons and neutrinos. The 
conservation of energy-momentum leads to the following 
set of time evolution equations for the photon and the 
baryon density contrasts S and velocity divergences 6 in 
the synchronous gauge [Zsj : 

= -^9^ - h, (2) 

Sb = -Ob - \k (3) 
for the density contrasts, and 

0^ = (^'^^'^ ^ '^■y^ + ane(JT{9b ~ ^7), (4) 

k = --Ob + cieS, + ^an,cTT{9^ - 9b), (5) 
a Spb 

for the velocity divergences. Here and throughout the 
paper, ax is the Thomson scattering cross section, 
is the electron number density, a is the scale factor, p 
is the background density, Cs is the sound speed given 
by Cs = 1/^/3(1 + R), ^ is the baryon-to-photon 
density ratio, cr^ is the photon shear, and the dot refers 
to the derivative with respect to the conformal time r. 
The variable h is the metric field in synchronous gauge, 
which evolves according to [i^ 

'^+^'' = -^{1) Pcr^%<5,(I + 3c,/), (6) 

where j S {f, 7, b, c} labels the different species, pcr is 
the critical density of the universe and flj = pj/pcr is 
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the ratio of the density of the j*'' species to the critical 
density. 

The tight-coupUng approximation allows us to set 6^ = 
Ob = Q-ybi with the photon-baryon velocity evolving as 



(1 + R%b = -RO^b + k^i^S-y - CT-,) 



clk^RSb, (7) 



and the photon density contrast evolving as [43 



5. + y:^s. 



Pels, = -| 



R 



1 + R 



h + h 



(8) 



Here, we have neglected the photon shear (tight-coupling 
regime) and the pressure term in Sb as it remains smaller 
than the term in 6-y prior to decoupling. Equation 
represents a driven harmonic oscillator with the compe- 
tition between gravitational infall and photon pressure 
giving rise to acoustic waves propagating in the photon- 
baryon fluid at the speed of sound. 

For the associated homogeneous equation, we look for 
solutions of the form 6-y oc expi Ludr' where w(t) is 
some phase function. The two solutions to the homo- 
geneous equation are simply sinfcr-s and cosfcrg, where 
rs(r) = /g Csdr' , the phase function is w = fccs, and we 
have made use of the WKB approximation. On large 
scales, the WKB approximation breaks down, but these 
modes are irrelevant for the BAO treatment as they only 
enter the horizon well after decoupling. The particular 
solution is constructed by integrating the driving term 
weighted by the Green's function of the two homogeneous 
solutions [43|. Thus, the time evolution of the acoustic 
waves in the photon component for all initial conditions 
prior to decoupling is given by 

(1 + i?)^/^(5-y(fc,T) = As sin fcrs(T) + Ac cos /cr^ (r) 

^ + R{T')f'^s\n [fcr,(T) - krs{T')]F{T')dT' , 

(9) 



-s Jo 



where As and Ac are determined by the initial conditions 
as described in [231, and 



where fg = an^aT is the differential optical depth. 
The Silk damping turns out to significantly affect both 
the shape and peak location of the BAO as we shall 
discuss later. After decoupling, the photons free stream, 
while baryons fall into the CDM potential wells under 
gravitational instability. Here, we only consider baryon 
fluctuations with wavelength larger than the Jeans scale. 

The above description of the density contrast evolution 
in fc-space can be intuitively and simply understood by 
looking at the evolution of the mass profile in the config- 
uration space [13, . The radial mass profile Mj of a 
species j, given by 



Mj{r,z)= Tj{k,z)^^^k^r^dk, (11) 
Jo '^^ 



where Tj{k,z) 
of the 



'5,ik,zf^dk, (12) 



5j{k,z)/k'^ is the transfer function 
J"" species, describes the redshift evolution of 
a point-like overdensity initially located at the origin. 
The location of the mass profile peak gives the physical 
radius of the spherical shell of the overdensity for a 
given species. For numerical computations, a Gaussian 
overdensity of width cr"^ is used instead of a point- like 
overdensity. This is done by multiplying the integrand 
of equation (dU by e~''^'^^/'^. 

Hereafter, we study the time evolution of the baryon mass 
profile for each mode in turn. We start from the well 
studied adiabatic case then move onto the isocurvature 
modes, since this will provide physical intuition into the 
effect of the isocurvature modes on the BAO. 



AD mode 



F(r) 



R 



3 \1 + R 



h + h 



(10) 



is the gravitational driving term which evolves differently 
for different initial conditions. Equation Q gives the 
time evolution of the photon density contrast in the tight- 
coupling regime. In this regime, the baryon density con- 
trast is related to its photon counterpart by Sb = jS,- 
On small scales, a correction to the tight-coupling ap- 
proximation must be applied when the Silk damping be- 
comes important, as photons leak out of overdcnse re- 
gions, dragging baryons with them. This is done by mul- 
tiplying the solution above by e"*^ ^'^^ , where the photon 
diffusion scale A: is given by 



CD 



1 R^ + A{l + R)/5 

~e (l + i?)2 ' 



The adiabatic mode is characterized by the requirement 
that the densities of all species are perturbed in propor- 
tion at some initial time such that 



(13) 



where the subscript i labels the initial time. Or equiv- 
alently, using the relative entropy between two species 
X and y given by S^y = - where Wx and 

Wy are the equation of state parameters of the species x 
and y respectively, we have that Sxy = for all pairs of 
species at the initial time. In addition, all velocity di- 
vergences are initially unperturbed. Therefore, using the 
initial conditions for the adiabatic mode [1^ , the photon 
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and baryon density contrasts are respectively given by 



V3 . 



k 



k /kj-j 



X / {1 + R{T')f^^ sm[kr,{T) ~ kr,{T')] 
Jo 

X F'^^(T')dr', (14) 
A k 

X [ {l + R{T')y/^sm[krs{T)^kr,{T')] 
Jo 

X F^^{T')dT'. (15) 

Thus, the adiabatic mode is only sourced by the grav- 
itational driving term F^^ . This driving term can be 
approximated by 

i^^^(fc,T)«2fc2c2jo(fcr,) 

on small and intermediate scales which reduces to 2A;^c^ 
at early times. On very large scales the above approx- 
imation breaks down, however, this does not affect our 
physical description of the BAO peak as these large-scale 
modes are well outside the horizon at decoupling and do 
not substantially influence the BAO features. The lack of 
an exact analytic expression for the driving term makes 
it difficult to derive exact analytic solutions for the time 
evolution of the photon and baryon density contrasts. 
Nevertheless, good approximations for the photon and 
baryon density contrasts are given by 



'2krsji{krs) x e 



(16) 



Therefore, at early times (fcrs(r) ^ 1) the density 
contrasts for the adiabatic mode, oc Si, oc {I — cos fcr,) 
couple to a cos kvg harmonic [1^ . 

Now, in the perfect tight-coupling approximation, that 
is if we omit the Silk damping correction in the density 
contrast equations, the baryon mass profile is given by 



Mb(r) cx (1 - H(r - r,)) ttT 



for r < rg. 
for r > rg, 



(17) 



which is obtained by substituting the density contrast ex- 
pression, without the Silk damping term, in the baryon 
mass profile expression. Here H{x) is the Heaviside step 
function. We observe that in the absence of Silk damp- 
ing, the baryon mass profile is quadratic at lower r and 
sharply peaked at a distance r(z) = ^^(z). This is illus- 
trated in Figure [T] where we show the effect of Silk damp- 
ing on the baryon mass profile at decoupling. We see that 
when we include the Silk damping term, the BAO peak 
is smoothed, attenuated and shifted to lower r. We can 
understand these features as follows. As we approach 
decoupling, the coupling between photons and baryons 




r(Mpc) 



FIG. 1: Effect of Silk damping on the baryon mass profile 
for the AD mode at decoupling. The solid curve represents 
the baryon mass profile without the Silk damping correction, 
while the dashed curve represents the baryon mass profile 
with the Silk damping factor turned on. In the absence of the 
damping term, the peak is located at r = = 144.5 Mpc. 



weakens and the photon mean free path increases. Pho- 
tons diffuse from overdensities to underdensities carrying 
baryons with them. Therefore, baryons leak out of the 
overdensity to both smaller and larger r, thereby smooth- 
ing and lowering the BAO peak. Due to the shape of the 
undamped mass profile (with no baryons on scales larger 
than the sound horizon). Silk damping has the effect of 
moving more baryons to larger scales. As a result, the 
BAO peak is at a slightly smaller distance than the sound 
horizon. As we will see later. Silk damping changes the 
shape of the mass profile for the adiabatic and isocurva- 
turc modes in different ways, due to the differing shapes 
of the undamped mass profiles. This has important con- 
sequences for our ability to use the BAO peak as a stan- 
dard ruler. 

The redshift evolution of the mass profile for the AD 
mode has previously been studied in the literature js^ ]. 
Initially the overdensities of all species coincide. As time 
evolves, the photon pressure drives acoustic waves in the 
photon-baryon fluid, while neutrinos free stream at the 
speed of light and the CDM remains at its initial lo- 
cation. In Figure [2] we show the redshift evolution of 
the CDM. baryon, photon and neutrino mass profiles. 
Prior to decoupling, photons drag baryons at the sound 
speed, leaving behind a void of baryons. Thus, the initial 
baryon point-like overdensity evolves in a spherical shell 
while the CDM overdensity collapses at the origin under 
gravitational instability, and the neutrinos free stream. 
After decoupling, photons free stream while baryons, free 
from the photons, collapse into the CDM potential wells. 
The baryon overdensity continues to collapse, pulling 
matter from the surrounding underdense regions to the 
overdense regions. As the baryon velocity divergence 
does not decay instantaneously at decoupling [SJ], the 
baryons only stall later at z ~ 500 with the consequence 
that the BAO peak is closer to 150 Mpc than 140 Mpc, 
the sound horizon size at decoupling. At z = 0, the 
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FIG. 2: Mass profile snapshots for the AD mode at difi'erent 
redshifts. The red, green, blue and purple curves respectively 
represent the CDM, the baryon, the photon and the neutrino 
mass profiles, (a) Well before decoupling (z = 3000), (b) At 
decoupling {z = 1080), (c) At late times {z — 0). The units 
of the mass profile are arbitrary but correctly scaled between 
panels. 



baryon mass profile displays two peaks, one near the ori- 



gin and a second peak at approximately 150 Mpc. 



B. NID mode 

The NID mode arises when the densities of the matter 
components are initially unperturbed while the initial 
perturbation in the neutrino density is balanced by its 
photon counterpart, keeping the curvature unperturbed. 
The initial perturbations are as follows: 



0, 6^ 



(18) 



These initial conditions imply that As = 0, thus excit- 
ing the cos kvs harmonic. The gravitational driving term 
contribution for this mode can be neglected without loss 
of accuracy, as the gravitational potential (related to h) , 
is initially unperturbed and only grows inside the hori- 
zon. This can also be understood by considering the 
right-hand side of equation ([6]). In the radiation domi- 
nated era, the photon and the neutrino density contrasts 
roughly cancel while the baryon and the CDM density 
contrasts remain small until the matter dominated era 
when they grow. The time evolution of the photon and 
the baryon density contrasts for the NID mode are given 

by 



•rNID 



\/3cs cosfcrs X e ^'^/^d ^ 



(19) 



cNID 



where Ri, 



(l - VScsCoskn) X 6-^^'=-, (20) 

Q^/flrad a-nd Rj — Cl~^/flrad s-^e respectively 
the fractional energy densities of neutrinos and photons 
at early times. The pressure due to an initial localized 
photon overdensity creates a baryon underdensity that 
propagates due to its coupling to photons and perturbs, 
through gravitational interaction, the CDM (see Figure 
[3]). In addition, isocurvature perturbations grow once 
they enter the horizon. It follows that the BAO peak in 
the case of the NID mode has smaller amplitude than in 
the adiabatic case. With time, the baryon and the CDM 
overdensities grow by pulling more matter from their 
surroundings, thus creating underdense regions around 
them. Note that the mass profile of a given species can 
be negative since the species can be initially perturbed 
positively, corresponding to an overdensity or negatively, 
corresponding to an underdensity, with respect to the 
background level. The final baryon mass profile displays 
a deeper trough between the two peaks compared to the 
adiabatic case. Most importantly, though the baryon 
overdensity in the NID mode evolves at earlier times like 
the baryon overdensity in the adiabatic mode as they 
both excite cos fcr^ harmonics, the final locations of the 
NID and the AD BAO peaks differ. At late times, the 
adiabatic mode becomes a superposition of sine and co- 
sine waves, departing from the NID mode and with the 
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FIG. 4: Effect of Silk damping on the baryon mass profile for 
the NID mode at decoupling. In the absence of the damping 
term (solid curve), the peak is located at r = = 144.5 Mpc. 
The dashed curve takes into account the effect of Silk damp- 
ing. 



damping. 

Figure S] shows the effect of Silk damping on the baryon 
mass profile at decoupling for the NID mode. In the 
absence of Silk damping, the AD and the NID BAO peak 
locations would coincide. The undamped baryon mass 
profile for the NID mode is given by 



Mb(r) cx (1 - H(r - r,.))r cx 



T for r < 
for r > Vg 



(21) 



As for the AD case, equation ((2T)) is obtained by omit- 
ting the damping factor in equation (|20p and substituting 
into equation (|12p for the mass profile. The baryon mass 
profile for the NID mode differs from the AD mode as 
it grows linearly with r until r(z) = rs(z) then falls to 
zero. For this reason, the shift in the BAO peak location 
due to Silk damping is larger than in the case of the AD 
mode for which, as previously mentioned, the undamped 
mass profile is quadratic in r for r < r^. The difference 
in the shape of the undamped mass profile also sets the 
difference in the width of the BAO peak. 



NIV mode 
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FIG. 3: Mass profile snapshots for the NID mode at different 
redshifts. The red, green, blue and purple curves respectively 
represent the CDM, the baryon, the photon and the neutrino 
mass profiles, (a) Well before decoupling {z = 3000), (b) At 
decoupling [z = 1080), (c) At late times {z = Q). The units 
of the mass profile are arbitrary but correctly scaled between 
panels. 



Unlike the other isocurvature modes, the NIV mode, like 
the AD mode, shows no relative entropy perturbation 
in the density field at some initial time. All the density 
perturbations are zero initially. The main difference with 
the AD mode is in the velocity field where the neutrino 
velocity divergence starts perturbed, being compensated 
by the photon-baryon velocity. The initial perturbations 
arc given by: 

R 

dc,i = 0, 9b^i — 9j^i — ~—9i,,i- (22) 



The NIV mode excites the sinfcr^ harmonic, so that 
undamped profile being convolved differently with Silk we can set Ac = in equation As in the case of the 
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FIG. 5: Mass profile snapshots for the NIV mode at difi'erent 
redshifts. The red, green, blue and purple curves respectively 
represent the CDM, the baryon, the photon and the neutrino 
mass profiles, (a) Well before decoupling (z = 3000), (b) At 
decoupling {z = 1080), (c) At late times {z — 0). The units 
of the mass profile are arbitrary but correctly scaled between 
panels. 



NID mode, the gravitational driving term contribution 



remains irrelevant at all times as all the densities start 
unperturbed and the perturbations only grow in the mat- 
ter dominated era. The time evolution of the photon and 
baryon density contrasts for the NIV mode are given by 



cNIV _ ^ cNIV 



sin fcr,,(T) X 6"'^'/'=^. (23) 



The non-zero initial velocity divergence of baryons and 
photons pushes the baryons and photons from the origin, 
thus creating an overdensity at approximately the scale 
of the sound horizon and a plateau at larger scales, 
in the baryon and photon mass profiles. The redshift 
evolution of the baryon overdensity for the NIV mode 
is shovifn in Figure [3 This is similar to the NID case, 
except that the baryon mass profile remains positive at 
all times due to the initial plateau. The BAO peak is 
at a different location as the sine harmonic convolves 
differently with Silk damping, compared to the cosine 
harmonic. 

In the absence of the Silk damping correction, the un- 
damped baryon mass profile for the NIV mode is given 

by 



M,(r) oc-^ In ^^-^ 
4 (r -f rsY 



cx 



for r <C r^, 
for r ^ 



r 



(24) 



The derivation of equation ([M)) is similar to the AD 
and NID cases. The undamped NIV mass profile grows 
quadratically with r for r < and peaks at r = as 
for the AD case. However, the shift in the BAO peak 
location due to Silk damping is not as significant as it is 
for the AD and the NID cases for the simple fact that the 
undamped mass profile does not abruptly fall off to zero 
after the peak as in the previous cases but decreases as 
before reaching a plateau of height proportional to 
Ts . This is due to the fact that the non-zero initial veloc- 
ity of photons carries baryons beyond the sound horizon. 
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FIG. 6: Effect of Silk damping on the baryon mass profile for 
the NIV mode at decoupling. In the absence of the damping 
term (solid curve), the peak is located at r = rs = 144.5 Mpc. 
The dashed curve takes into account Silk damping. 
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compared to if they started from rest. Figure [6] shows 
the effect of Silk damping on the undamped baryon mass 
profile. In contrast to the AD and NID cases, the BAO 
peak is slightly shifted to higher r. 



D. CI &: BI modes 



The CI and the BI modes have been well studied in the 
literature [ill, HI] . The CI and BI modes arc similar 
in that the perturbation starts in the CDM density con- 
trast and the baryon density contrast respectively while 
the other species are initially unperturbed. This can be 
written at some initial time as 



Set = 1, Sb 



for the CI mode, and as 



0, 



(25) 



Sc,i = 0, 6b,i — 1, 6j^i — 5u^i — 0, 



(26) 



for the BI mode. The CI and BI initial conditions dic- 
tate that As = ~^^^c,Q for the CI mode and As = 
~^^^b,Q for the BI mode, while Ac = in both cases, 
thus exciting the sinfcrg harmonic [1^ H^. The con- 
stants f2c,o and 57f,,o are respectively the CDM and the 
baryon densities today. The driving term is negligible 
in the radiation domination era as the photon and the 
neutrino densities are initially unperturbed but becomes 
important in the matter domination era as the matter 
perturbation sources the gravitational potential [42 1. 



The time evolution of the photon and baryon density 
contrasts for the CI and BI modes is given by [i^l 



xCi 



3""^"' k 

V3 r 



i^c,o "7^ sin A;rs(T) X e ^^^^^^ 



^ ^ (1 + R{r')f'^ sin [fcr,(r) - fcr,(r')] 



"b 



90 ^ 

3^3 r 



(27) 



a/3 

-2ric.on~ sin krs{T) 
k 



4 k 

X F^'^(r')rfr' X e 



{1 + R{t')Y'^ sin [kr,{T) - fcr,(r')] 

(28) 
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FIG. 7: Mass profile snapshots for the CI mode at different 
redshifts. The red, green, blue and purple curves respectively 
represent the CDM, the baryon, the photon and the neutrino 
mass profiles, (a) Well before decoupling [z = 3000), (b) At 
decoupling (z = 1080), (c) At late times (z = 0). The units 
of the mass profile are arbitrary but correctly scaled between 
panels. 
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for the CI mode, and by 



+ ^ /V + i?(T'))'/' sin [fcr,(r) - fcr,(T')] 



3^3 r 



(29) 



\/3 

1 - 2f7h,o^ sin fcrs(T) 



4 /fc 



(1 + i?(T'))^/^ sin [A:r,(T) - fcr,(r')] 



(30) 



for the BI mode. Equations (|27ll30p are exact but require 
a perfect knowledge of the gravitational driving term. 
This makes the derivation of simple explicit analytic ex- 
pressions for the CI and BI modes harder as compared 
to the AD, NID and NIV modes. Therefore, we do not 
discuss the effect of Silk damping on the BAO peak fea- 
ture for these modes. However, one thing to notice is the 
k^^ dependence of the baryon density contrast for the 
CI and BI modes that washes out perturbations on small 
scales while amplifying them on large scales. This redis- 
tribution of power results in a flattening of the baryon 
mass profile for these modes. On small scales, the 
and the Silk damping factors have similar effects on the 
BAO peak as they both suppress perturbations on these 
scales. However there are two main differences. Firstly, 
Silk damping does not act on large scales while the k~^ 
factor amplifies large scale perturbations. Secondly, Silk 
damping only becomes significant around recombination 
while the k^^ factor redistributes the power at all times, 
hindering the development of a well defined BAO peak 
but producing a knee instead. 

For the CI mode, an overdensity in the CDM compo- 
nent tends to affect, through gravitational attraction, the 
baryon density component by gathering baryons into an 
overdensity but the photon pressure opposes this process 
until decoupling. 

One should note that an initial overdensity in the 
photon component would more easily affect the baryon 
component than an initial overdensity in the CDM 
component, the reason being the high photon pressure 
at earlier times. Therefore, the perturbation takes longer 
to imprint ripples onto the homogeneous sea of baryons. 
Figure [7] represents the time evolution of the baryon 
mass profile for the CI mode. Prior to decoupling, 
the CDM overdensity grows but does not significantly 
affect the baryon component. After decoupling a baryon 
overdensity develops through gravitational interaction 
with the CDM but fails to display a well defined BAO 
peak. 

For the BI mode, an initial overdensity in the baryon 
component affects the CDM component through gravi- 
tational attraction, but does not significantly grow due 
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FIG. 8: Mass profile snapshots for the BI mode at different 
redshifts. The red, green, blue and purple curves respectively 
represent the CDM, the baryon, the photon and the neutrino 
mass profiles, (a) Well before decoupling {z = 3000), (b) At 
decoupling {z — 1080), (c) At late times {z — 0). The units 
of the mass profile are arbitrary but correctly scaled between 
panels. 



to the photon pressure at earlier times that tends to 
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widen and even wash out the baryon overdensity as 
can be seen in Figure [S] With a similar process as for 
the CI mode, the overdensity becomes a knee at late 
times. Although the Silk damping affects the CI and BI 
modes, its effect is not as si gnif icant as in previous cases 
(for a discussion of this see |46l . |48| ) . We recall that the 
Silk damping tends to suppress power on small scales 
while these modes are already significantly reduced by 
the factor for the CI and BI modes. In addition, 
the fact that the CI and BI modes fail to display a well 
defined BAO peak makes less noticeable the effect of 
Silk damping on the BAO peak. 



sitions. After decoupling, though the BAO peak position 
still increases until z sa 500 due to the bulk velocity, the 
separation between the modes remains constant until to- 
day. 

III. IMPACT OF ISOCURVATURE MODES ON 
DARK ENERGY CONSTRAINTS 

The aim of this section is to quantify the potential im- 
pact of isocurvature modes on dark energy studies based 
on current- & next-generation datasets. 



Statistical Formalism 



E. Time evolution of the BAO peak position 

We saw in previous subsections that in the absence of 
Silk damping, the BAO peak location for all the modes 
would coincide at all times as the acoustic wave in the 
photon-baryon fiuid propagates at the same sound speed 
irrespective of the initial conditions. Here we consider 
the effect of Silk damping on the evolution of the BAO 
peak location for different modes. 




1000 1500 
Conformal time t 



FIG. 9: Time-evolution of the baryon mass profile peak loca- 
tion (see Figures [1 [3] and [51 for the AD, NID and NIV modes. 
These curves were obtained numerically from the evolution of 
the mass profile curves. The dashed vertical line indicates the 
epoch of recombination. 

Figure IH] shows the time evolution of the BAO peak 
position for the AD, NID and NIV modes. We do not 
include the CI and BI modes as they fail to display a 
defined BAO peak. At early times, the BAO peak po- 
sitions for the AD, NID and NIV coincide as the Silk 
damping factor equals one at early times. With time, the 
damping scale fe^^ increases as the photon-baryon cou- 
pling weakens, the three modes depart from each other 
and the separation increases up to decoupling. At de- 
coupling, ku ss 0.15 Mpc~^, leading to a separation of 
about 15 Mpc between the NID and NIV BAO peak po- 



A convenient way of quantifying the accuracy with 
which cosmological parameters can be measured from a 
given dataset is the Fisher matrix formalism (see [ll| 
for a review). If x is our observable (the CMB or the 
galaxy power spectrum in our case), it can be modeled 
as a N-dimensional random variable whose probability 
distribution /(x;0) depends on a vector of cosmological 
parameters 9 that we wish to estimate, ^x; 0) is also 
known as the likelihood of observing a set of data given 
a model characterized by 6. 

Ill this study, we consider a spatially flat cosmological 
model described by the following parameters: the baryon 
density Wb, the CDM density ujc, the density of the dark 
energy component fix, the optical depth r, the spectral 
index rig and the scalar amplitude As . We allow for dark 
energy models that vary with time and parametrize the 
dark energy equation of state as w{a) ^ wq + {I — a)wa 
jsvl [ssj where a = 1/(1-1-2;) and wq and Wa are included 
in the parameter space. For the isocurvature modes, we 
adopt the parametrization implemented in [39j . where 
the AD, CI, NID, NIV modes and their cross-correlations 
are described by 10 parameters, Zij, measuring the frac- 
tional contributions of the various correlations (auto and 
cross) to the total power spectrum. We do not consider 
the BI mode as it has the same spectra as the CI mode. 
In terms of these fractional parameters, the total isocur- 
vature fraction /igo is given by 



/iso 



^ISO 



ZlSO + 2(AD,AD) 



(31) 



1 



-(AD, AD) 



is the total isocurvature 



where zjso = 
contribution. 

Defining the auto- and cross-correlated primordial 
power spectra as follows 

P„-(fc) = (32) 

and the observed full power spectrum that includes 
adiabatic and all different modes as 



P{k) = J2 A^,lP^Ak), 



(33) 
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the spectral indices of the cross correlated modes are 
given by = with their amplitudes Aij oc z^ . 

The constraint X^l "=1 ^fj = 1 requires that the isocur- 
vature parameters Zij exist on the surface of a 9 di- 
mensional sphere of unit radius. For further details on 
how this parametrization relates to others in the litera- 
ture, see [22|. The parameter As rescales the unit power 
CMB temperature spectrum to its usual amplitude as, 
Ci ^ ISOOO/iK^ where Ce is the fiducial CMB 

temperature spectrum with unit power. 

Fiducial model 



As 



Wo Wa 



0.02205 0.12495 0.7 



0.1 1.0 15.7 -1.0 0.0 



TABLE I: Values for the parameters of the fiducial cosmolog- 
ical model. 



The Fisher matrix is defined by 



92 ln;(x; 9) 



(34) 



The Cramer- Rao inequality shows that F,"^ is the small- 
est variance possible for an unbiased estimator of the 
parameter 9i. In this case, is the most optimistic 
covariance matrix of the dataset [llj and the forecasted 
error bar for 6i is 



[F- 



(35) 



The fiducial model 6 around which the Fisher matrix is 
computed is chosen to be a ACDM universe with adia- 
batic initial conditions. The cosmological parameter val- 
ues for the fiducial model are given in table H] We use the 
CMB forecasts from the Planck experiment in addition 
to each LSS data set and compute the full Fisher matrix 



pCMB 
ij 



pLSS 



for the cosmological parameter set. 



1. Large Scale Structure (LSS) surveys 

Over the next decade, the increase in the number and 
quality of data from LSS surveys will drive fundamental 
improvements in precision cosmology. As these galaxy 
surveys cover increasingly larger volumes, they will pro- 
vide unprecedented probes of scales at which significant 
cosmological information is available. 

The potential of the BAO method as a powerful source 
of cosmological information has been recognized and 
measuring the BAO peak at multiple redshifts is now 
regarded as the primary science of major future LSS 
surveys. We consider two such BAO experiments, one 
of which is the Baryon Oscillation Spectroscopy Survey 
(Boss). Boss will measure the redshifts of 1.5 million 
luminous red galaxies (LRGs) over a quarter of the sky 



to a depth of z = 0.7. In addition to being a redshift sur- 
vey. Boss will be the first attempt to resolve the BAO 
peak in the high-z density field (2 < z < 3), as allowed 
by mapping absorption lines from neutral hydrogen, in 
the spectra of 160 000 distant quasars [20j . 

The Advanced Dark Energy Physics Telescope 
(Adept) is the second future LSS survey that we con- 
sider. It is a space-based experiment aiming at mapping 
galaxies in the redshift range 1 < z < 2 and over 28, 600 
sq. deg. of the sky [Toj . 

The BAO peak manifests as oscillations in the matter 
power spectrum with the size of the sound horizon deter- 
mining the frequency of these oscillations. However, the 
matter power spectrum is a rich statistic whose features 
at different scales provide specific cosmological informa- 
tion. The matter power spectrum is defined as 

P(fc, z) = D{z)^Pp„rnik)T\k) (36) 

where D{z) is the growth rate of structure, Pprim{k) is 
the primordial power spectrum and T{k) is the transfer 
function. The first source of information is the baryon 
acoustic oscillations, with their wavenumber k ~ 27r/rs 
being set by the size of the sound horizon at decou- 
pling r^. Since this characteristic scale is calibrated 
by the CMB, measuring the wavelength of these oscil- 
lations both in the radial and tangential directions deliv- 
ers Da{z) and H{z) respectively. The overall shape of 
the matter power spectrum is a second source of infor- 
mation. Any features which deviate from a power law, 
such as the turnover, provides an additional characteris- 
tic scale which is required by the Alcock-Paczynski test 
to be isotropic @. Lastly, the overall time evolution of 
the amplitude informs us about D(z), the growth rate of 
structure. 

In reality we measure the power spectrum as mapped 
by galaxies which are biased tracers of the underlying 
matter distribution. We can write the galaxy power spec- 
trum as Pg{k,z) = b{k, z)'^P{k, z) where b{z,k) repre- 
sents this bias resulting from the effects of galaxy forma- 
tion and evolution. On the scales of the BAO, the bias 
can be regarded as smooth, i.e., b{z, k) = b{z). Any scale 
dependence that is not taken into account is not likely to 
lead to oscillations in Fourier space (TJ]. Furthermore, 
the galaxy power spectrum measured in redshift space 
is distorted relative to the power spectrum in real space 
as a result of galaxy peculiar velocities. Because galax- 
ies moving towards an overdensity along the line of sight 
appear further away than equidistant galaxies moving 
in the tangential direction, structures appear " squeezed" 
in redshift space, with the amount of the distortion de- 
termined by the growth rate. On large scales this has 
been shown to give rise to an angle-dependent distortion 
which leads to a multiplicative change in the power that 
is a function of angle, ie. Pg^p = (l + f3{z)fi'^) Pg{k,z) 
where ji is the angle with respect to the line of sight and 
/3 = f /b where 



/ 



ainL>(a) 
dlna 



0.6 



(37) 
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We follow [l5| and include a free shot noise parameter as 
well as a redshift distortion parameter /3(z) in each bin. 

Assuming the likelihood function of the band powers 
of the galaxy power spectrum to be Gaussian, the Fisher 
matrix can be approximated as [ill, '■ 



kmin 
1 pkrnax 



-1 Jk 



dp, dpj - ' 2{27Ty 

dill P{k,n) dlnP{k,fi) 



dpi 
2{2Tif 



dpj 



where, 



Veff{k,n) 



ngPg{k){l+l3ii 



2\2 



V, 



(38) 



(39) 



Boss 




z 


kmax/ii MpC ^ 


b 


V/ Gpc^ 


Area/deg^ 


3 X 10"^ 


z < 0.35 


0.09 


2.13 


0.74 


10,000 


3 X lO"-* 


0.35 < 2 < 0.6 0.11 


1.25 


2.83 


10,000 


3.1 X 10^2 


2 < 2 < 3 


0.3 


-0.18 


2.48 


6000 


Adept 




z 




b 


V/ Gpc^* 


Area/deg^ 


3 X 10"'' 


Kz < 1.25 


0.16 


2.97 


17.7 


28600 


3 X lO"-* 


1.25 < z < 1.5 


0.11 


3.21 


19.7 


28600 


3 X 10"^ 


1.5 < z < 1.75 


0.22 


3.44 


21.0 


28600 


3 X 10"'' 


1.75 <z <2 


0.25 


3.67 


21.7 


28600 



TABLE II: Table summarizing the survey parameters for 
Boss and Adept, for different redshift bins (centered at the 
middle of the redshift bin). 



f is the unit vector along the line of sight and k is the 
wave vector with norm k = \k\. Here V is the survey 
volume contained in a given redshift bin and ng(f) is the 
selection function of the survey, dictating the a priori 
expectation value for the comoving number density of 
galaxies. We take this to be a constant. Vefj is the 
effective volume of the survey and takes into account the 
impact of the shot noise from undersampled regions [l^l ■ 

The derivatives of the power spectrum with respect to 
the cosmological parameters in table U and to the isocur- 
vature parameters are respectively shown in Figures 1161 
and [17] in the appendix. 

The Boss and Adept survey parameters are summa- 
rized in Table ini Note that the value kmin is always taken 
as the lowest possible and has been shown to have a neg- 
ligible effect on the error forecasts. The smallest scale 
included, given by kmnj^ in the analysis does however 
impact on the results [la, [13 • Following [l^ we adopt 
conservative values for kmax by requiring a{R) = 0.2 at 
a corresponding R = where a{R) is defined similarly 
to the normalization erg = a{R = 8/i~^Mpc), but for 
a general scale R. For the BOSS z = 3 bin, we restrict 
kmax = 0.3 h Mpc~^. Although the application of the 
Fisher matrix formalism to the case of a Lyman-a forest 
is different to that of a galaxy survey in various subtle 
ways, the measurement of the BAO scale from a Lyman- 
a survey has however been explored extensively in [isj . 
Assuming typical values, we follow their prescription to 
estimate the shot noise contribution (l/ng)for the z = 3 
bin of the Boss survey. Note that b for the z = 3 bin is 
our estimate the linear bias between the flux and density 
field. 



2. Cosmic microwave background (CMB) surveys 



Non-adiabatic initial conditions lead to very distinct fea- 
tures in the temperature anisotropics, with isocurvature 
modes producing acoustic oscillations that are out of 
phase with the adiabatic mode and hence a set of peaks 
in the temperature anisotropy power spectrum that are 
slightly shifted. Furthermore, CMB polarization pro- 
vides a robust signature of isocurvature perturbations 
[2i |. The latest WMAP data has confirmed that the 
initial perturbations were mainly of adiabatic type (soj 
with the possible presence of a subdominant isocurva- 
ture contribution, which could be detected in future high- 
precision experiments such as Planck 0. The higher 
resolution of Planck over WMAP will allow for the mea- 
surement of the CMB power spectrum on much smaller 
scales and the use of 9 observational bands will improve 
the modeling of astrophysical foregrounds. 

We follow the analysis in (Z^l and model the Planck 
dataset as CMB temperature and polarization maps of 
80% of the sky measured in the two frequency bands 
where the CMB signal dominates. The details of the 
experiment are given in Table Hill The maps are taken to 
have no foreground contribution, assuming that the other 
frequency channels can be used to remove them. The 
remaining 20% of the sky is assumed to be contaminated 
by galactic emission. We exclude polarization data at 
£ < 30 in order to weaken the forecasted constraint on the 
optical depth to cr(T) = 0.01 in agreement with studies 
that include foreground modeling p?5l |. 



v/ GHz At(mK) Ap(mK) 



2000 2500 143 
217 



8' 5.2 
5.5' 11.7 



10.8 
24.3 



The CMB data primarily provides information about 
the initial conditions of our Universe in this analysis. 



TABLE III: Summary of the experiment specifications for 
Planck. 
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For the CMB, the Fisher matrix is computed using 

(40) 



CMB 



\ ^ \ ^ dCxf. 1-1 dCxi 



X,Y 



dpj 



where Cxi is the power in the 1 muhipole for X 
T, E, B given by 



[Cov^; 



XX — 



{Cxi + f^l) 



(41) 



where N(, the noise level, depends on the data type. The 
noise is specified by the experiment. 

Because there is a strict geometric degeneracy be- 
tween r^A, li'o and Wa, finding the derivatives of the 
dark energy equation of state (EOS) parameters while 
keeping f^A fixed artificially breaks this degeneracy. To 
this end, we follow [43| and start with computing the 
Fisher matrix for the CMB with the following parame- 
ters: p^jojf,, ojc, ^s, T, fis, where 9s is the angular size 
of the sound horizon. This can be written as 



_ ^ s (^cmfe) 
^{^cmb ) 



(42) 



where rs{zcmb) is the sound horizon given in equation [T] 
and r{zcmb) is the comoving distance to the last scatter- 
ing surface 



r{z) 



H{z') 



-dz'. 



(43) 



To compute the derivative dCi /dOg , we use the transfor- 
mation 



dCe 



(44) 



and when evaluating -^rf-, ujc and iOh must stay fixed by 



compensating with h through 



(45) 



The resulting Fisher matrix F is then trans- 
formed back into F, corresponding to the parameters 
p'={wb, ^^A, T, ris, As,wo,Wa} using 



F. 



E 

n.rn 



dprn 
dpi 



■Frr 



dpn 

'dp 3 



(46) 



The non-trivial expressions needed for the Jacobian are 
derivatives of 9s with respect to wq, Wa, i^A, and uJh- 

For the calculations of the derivatives of the power 
spectrum with respect to the isocurvature amplitudes in 
equation[Ml we have adopted the treatment in [2^ where 
the pure isocurvature modes are normalized to have the 
same power in their CMB temperature spectra as the 
adiabatic model. This normalization is applied to both 
the CMB and LSS spectra. 



Experiment 


Planck 
+ Boss 




Planck 

-1- Adept 




Parameters 


cr{wo) 


cr{wa) 


(j(u)o) 


O-(w'a) 


Adiabatic mode 


0.045 


0.083 


0.035 


0.060 


Adiabatic -I- 1 ISO mode 


ad-|-ci-|-(ad,ci) 


0.058 


0.087 


0.042 


0.063 


ad-|-nid-|-(ad,nid) 


0.048 


0.086 


0.038 


0.062 


ad+niv+(ad,niv) 


0.049 


0.091 


0.037 


0.070 


Adiabatic -|- 2 ISO modes 


AD-f Cl-f NlDH-corr 


0.058 


0.089 


0.043 


0.065 


AD-f Cl-fNlv-f corr 


0.069 


0.11 


0.055 


0.090 


AD+Nlv-l-NlD-l-corr 


0.055 


0.098 


0.043 


0.080 


Adiabatic -|- all ISO modes 0.073 


0.12 


0.061 


0.10 



TABLE IV: Table summarizing tiie constraints on (wojWa) 
for adiabatic and admixtures of uncorrelated adiabatic and 
isocurvature modes, marginalizing other all other parame- 
ters, for the Boss and Adept experiments. The decrease in 
the figure of merit when all isocurvature modes are added to 
the adiabatic case are FoM/so=0.54 FoMad and FoM/so = 
0.47 YoMad for Boss-I-Planck and ADEPT-f Planck respec- 
tively. The fiducial model assumes adiabaticity. 



B. The impact of isocurvature modes on dark 
energy 

In this section we consider the impact of admitting 
isocurvature initial conditions on the constraints on the 
dark energy parameters. We follow [i^l and choose not 
to focus on the constraints on Q.x- 

We compute the potential errors on wq and Wa for dif- 
ferent subsets of adiabatic and isocurvature initial con- 
ditions while marginalizing over all other cosmological 
parameters. The results for both the Boss and Adept 
experiments are summarized in table IIVI We find a sys- 
tematic degradation of the viable constraints on dark en- 
ergy as more degrees of freedom are added. In order to 
quantify the constraining power of the data, we compute 
the Dark Energy Task Force (DETF) Figure of merit 
(FoM), which is defined as the reciprocal of the area in 
the wq — Wg plane, enclosing the 95% confidence limit 
(CL) region [40|. We are concerned with the change in 
the FoM when isocurvature modes are introduced rela- 
tive to the case of pure adiabaticity. The BOSS FoM is 
found to decrease by 46% from pure adiabaticity to the 
case in which all isocurvature modes arc admitted in ad- 
dition to the adiabatic, while the Adept FoM degrades 
by 53%. We note that the results quoted here arc slightly 
different to those quoted in Q. This is due to the differ- 
ent normalization method used in Q , which follows [3^ , 
whereas we follow [2^. However, the different normal- 
ization methods used have little impact on the results 
reported here, which display a similar trend to @ but 
with a slightly larger relative degradation of parameter 
errors when isocurvature modes are included. 
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The results suggest that no single mode in particu- 
lar and its correlation are responsible for the change in 
the allowable (wQ^Wa) region, but rather a mixture of 
all extra degrees of freedom. In order to determine the 
combination of parameters that is responsible for this 
degradation, we diagonalize the full 17x17 Fisher matrix 
corresponding to the Planck and LSS datasets sepa- 
rately and find the eigenvector with the smallest eigen- 
value, corresponding to the direction that is least con- 
strained by the data, that is, we consider the worst case 
degenerate scenario. Considering the Planck data alone 
and discarding the degeneracy in the wq — Wa direction 
which is the main degenerate direction, four parameters, 
namely As, ^(ad,niv), ^(AD,nid) and ^(nid,nid) de- 
fine the most degenerate direction involving isocurvature 
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FIG. 10: Main contributions to the degenerate direction with 
the highest isocurvature fraction (a) in the CMB data from 
Planck alone, (b) in the matter power spectrum using the 
Boss dataset alone. The red solid line is the total derivative in 
the considered degenerate direction, which cancels to within 
the limits allowed by Planck and Boss error bars (yellow 
region). For the matter power spectrum, we only plot the 
yellow region up to kmax = 0.1 h Mpc~^ in this redshift bin. 
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FIG. 11: (a) The galaxy correlation functions for the Boss 
survey for different redshift bins. Sohd lines represent the 
purely adiabatic ACDM fiducial model while dashed lines 
represent a mixed model with /jgo = 95%, wq — —1.07 & 
Wa = —0.133. z = 0.35 (black), z = 0.6 (blue) and z — 3 
(red), (b) Relative difference of the fiducial and mixed mod- 
els. 



modes, with the impact of the scalar amplitude being 
compensated for by a combination of the cross-correlated 
modes and the optical depth. The degenerate direction 
in the LSS data (using BOSS as an example) is more 
complicated and involves a combination of isocurvature 
parameters and dark energy parameters, namely Qx, ^s, 

Wo, Wa, ^(AD,CI>j ^(AD,NIV>, ^(AD,NID> and ^(CIjNIV)- 

Figure [10] shows how the perturbations in the different 
parameters contribute to the the total change in the CMB 
power spectrum and the matter power spectrum. 

The total derivative (shown in red) lies within the 
noise limits of the respective experiments, making the 
net change undetectable by the data. 

Clearly the dark energy model is degenerate with the 
particular combination of isocurvature modes in the B AO 
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Experiment 


Planck + Boss 




Planck + Adept 




Parameters 




SWa 


Swq 


5Wa 


Adiabatic + 1 ISO mode 


ad+ci+(ad,ci) 


-0.12 (2.7) 


0.077 (0.9) 


-0.076 (2.2) 


0.072 (1.2) 


ad+nid+(ad,nid) 


-0.060 (1.3) 


0.084 (1.0) 


0.050 (1.4) 


-0.067 (1.1) 


ad+niv+(ad,niv) 


0.091 (2.0) 


-0.17 (2.0) 


-0.055 (1.6) 


0.16 (2.7) 


Adiabatic + 2 ISO modes 


AD+CI+NID+corr 


-0.11 (2.4) 


0.094 (1.1) 


-0.075 (2.1) 


0.073 (1.2) 


AD+Cl+Nlv+corr 


-0.19 (4.2) 


0.22 (2.7) 


0.16 (4.6) 


-0.23 (3.8) 


AD+NIV+NID+COrr 


-0.11 (2.4) 


0.19 (2.3) 


-0.088 (2.5) 


0.20 (3.3) 


Adiabatic + all ISO modes 0.22 (4.9) 


-0.29 (3.5) 


0.20 (5.7) 


-0.29 (4.8) 



TABLE V: Table summarizing the biases on (wojWa) that could arise from the incorrect assumption of adiabatic initial 
conditions, given a universe with an admixtures of uncorrelated adiabatic and isocurvature modes for the BOSS and Adept 
experiments. The quantities in brackets are the biases, quoted in number of la error bars corresponding to the case when pure 
adiabaticity is assumed. 



data. The implication is that the constraints on dark 
energy are at risk of being substantially biased if adia- 
baticity is incorrectly assumed. To emphasize this point, 



Figure 11(a) compares the correlation function, defined 

by 



kr 



(47) 



where P{k) is the matter power spectrum, that would 
be measured today for our fiducial ACDM model assum- 
ing pure adiabaticity, to a cosmological model assuming 
dynamical dark energy, described by wq = —1.07 and 
Wa = —0.133, and an admixture of initial conditions, 95% 
of which is isocurvature in nature. The correlation func- 
tion is degenerate in all three redshift bins of the BOSS 
experiment. Note that this degeneracy is completely bro- 
ken by the CMB data. We now wish to quantify this bias. 

For a Gaussian-distributed likelihood function, it can 
be shown that the linear bias in a set of parameters that 
we wish to constrain, 56 due to erroneous values of a 



set of fixed parameters. 



56, 



is \M 



L J im rnj 



(48) 



where F^^ is the Fisher sub-matrix for the parameters 
we wish to constrain and F^"^ is a Fisher sub-matrix 
constructed from the product of the derivatives of the 
power spectrum with respect to the parameters being 
constrained and those which arc being fixed. In our case 
i labels the isocurvature mode amplitudes, incorrectly 
fixed to zero, m labels the eight cosmological parameters 
that are biased, and i labels the subset of two dark energy 
parameters whose bias is of interest to us. In order to 
set 5(j)j, we diagonalize the combined Planck and large- 
scale structure (LSS) Fisher matrix and select the eigen- 
vector, with the smallest eigenvalue A^. This corre- 
sponds to the direction in parameter space which is least 



constrained by the data. We then take 5(1) j = ^ ^-^ji 

where M depends on the total number of cosmological 
and isocurvature parameters. 

Wc first consider the case of an admixture of the adi- 
abatic mode and the CDM isocurvature mode. For this 
case we find the biases in the dark energy parameters to 
be 5wo = -0.12 and 5wa = 0.077 for the Boss experi- 
ment. Comparing the mean biases to the la constraints 
obtained when pure adiabaticity is assumed, we find that 
neglecting this isocurvature contribution leads to a 2. 7 a 
and 0.9(7 error in the dark energy parameter estimates 
for wq and Wa respectively, when compared to the error 
forecasts assuming adiabaticity. If we repeat the calcula- 
tion for the more advanced experiment Adept, we find 
5wo = —0.076 and 5wa = 0.072, equivalent to 2.2a and 
1.2a errors in the dark energy parameters respectively 
when compared to the adiabatic constraints. 

Although no theoretical models for generating the neu- 
trino isocurvature models have thus far been proposed, 
we would like to conduct a comprehensive exploration 
of the impact of the initial conditions on the BAO con- 
straints and therefore admit all possible isocurvature de- 
grees of freedom. Table |V] summarizes the biases for dif- 
ferent admixtures of adiabatic and subsets of isocurva- 
ture modes. The results are consistent with the degener- 
ate directions in parameter space identified earlier. For 
example, the admittance of the CI and NIV isocurvature 
modes and their cross correlations has the potential to 
cause a bias in the dark energy parameters by as much 
as 4.2(7, in the case of BOSS. 

For the case of an admixture of adiabatic and all isocur- 
vature modes and their cross correlations, wc find that 
the biases are 5wo = 0.22 and dwa = —0.29 for the 
Boss experiment. This means that if the initial con- 
ditions of our universe are comprised of a sub-dominant 
contribution from all isocurvature modes (within the la 
constraints from the Planck and BOSS experiments), 
the assumption of adiabaticity could lead to an incorrect 
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4.9<T detection of non-A dark energy model or a 3.5a false 
claim of dynamics. Alternatively, A could be found to be 
consistent with the data when in fact w{z) 7^—1. The po- 
tential bias incurred by the adiabatic assumption in the 
case of the Adept experiment has a mean of Swq = 0.20 
(equivalent to b.7a) while the measurement of Wa could 
be inaccurate at the level of only 4.8ct. 



1. Constraints on isocurvature modes from the LSS data 

We now consider the impact of the large scale struc- 
ture information on isocurvature constraints. Although 
allowing for isocurvature modes degrades the dark en- 
ergy constraints relative to the pure adiabatic case, this 
analysis has revealed a powerful positive. The forecasted 
errors on the isocurvature parameters based on the CMB 
data alone and in conjunction with the LSS experiments 
are compared in tables IVIl IVIII and I Villi respectively for 
single, double and fully correlated isocurvature modes. 
We find that the error bars on the isocurvature parame- 
ters decrease by a few percent to as much as 96% for cer- 
tain modes when the LSS data (either BOSS or Adept) 
is added to the Planck data. Assuming an adiabatic 
fiducial model, the measurement of the BAO in the first 
rcdshift bin of the BOSS experiment and the CMB by 
Planck will reduce the allowed isocurvature fraction 
from 5.6% for the CMB data only to 3.5%, and to 3.1% 
and 3% when adding the information from LSS in the 
second and the third redshift bins. 



Adiabatic -I- 1 ISO mode 


Planck -i- Boss 


Planck + Adept 


(ad,ad) 0.31 (4) 


0.32 (1) 


{ci,ci) 0.012 (5) 


0.012 (6) 


(ad,ci) 0.035 (96) 


0.032 (96) 


(ad, ad) 0.027 (14) 


0.26 (18) 


{nid,nid) 0.0050 (6) 


0.0049 (9) 


(ad,nid> 0.016 (13) 


0.015 (18) 


(ad,ad) 0.25 (21) 


0.24 (24) 


(niv,niv) 0.0089 (5) 


0.0089 (6) 


{ad,niv) 0.090 (90) 


0.057 (94) 



TABLE VI: Forecasted uncertainties on isocurvature param- 
eters for diflerent cases for the Planck and LSS data (Boss 
and Adept) for single isocurvature modes. The percentage 
improvement in la errors wlien tlie LSS data is added to the 
Planck data is shown in brackets. 

The reason for this stems from the fact that the con- 
sidered degenerate direction in parameter space for the 
CMB data ((ad,niv), (ad,nid), (nid,nid). As) differs 
from the degenerate direction of the LSS data ((ad,niv), 

(AD,NID), nx, Wo, Wa, As). 

Figure [T^] illustrates the different directions by show- 
ing the la error ellipses for the main isocurvature con- 



Adiabatic -I- 2 ISO modes 




Planck + Boss 


Planck + Adept 


(ad, ad) 


0.36 (19) 


0.34 (24) 


(ci,ci) 


0.20 (5) 


0.020 (7) 


(nid.nid) 


0.0072 (2) 


0.0071 (3) 


(ad,ci) 


0.035 (96) 


0.032 (96) 


(ad,nid) 


0.029 (8) 


0.027 (13) 


(ci,nid) 


0.015 (2) 


0.015 (3) 


(ad, ad) 


0.43 (4) 


0.43 (4) 


(ci,ci) 


0.017 (5) 


0.017 (5) 


(niv,niv) 


0.0095 (4) 


0.0095 (4) 


(ad,ci) 


0.044 (96) 


0.043 (96) 


(ad,niv) 


0.096 (91) 


0.068 (94) 


(CI,NIV) 


0.017 (3) 


0.017 (3) 


(ad, ad) 


0.28 (38) 


0.27 (40) 


(nid,nid) 


0.0083 (3) 


0.0082 (4) 


(niv,niv) 


0.018 (3) 


0.018 (3) 


(ad,nid) 


0.028 (10) 


0.028 (11) 


(ad,niv) 


0.093 (90) 


0.062 (93) 


(nid,niv) 


0.016 (3) 


0.016 (4) 



TABLE VII: Forecasted uncertainties on isocurvature param- 
eters for different cases for the Planck and LSS data (Boss 
and Adept) for double isocurvature modes. The percentage 
improvement in la errors when the LSS data is added to the 
Planck data is shown in brackets. 



Adiabatic + all ISO modes 

Planck + Boss Planck + Adept 



(ad, ad) 


0.51 (13) 


0.51 (14) 


{CI,Cl) 


0.047 (18) 


0.044 (23) 


(nid,nid) 


0.017 (11) 


0.016(14) 


{niv,niv) 


0.039 (17) 


0.037 (22) 


(ad,ci) 


0.061(94) 


0.057 (95) 


(ad,nid) 


0.088 (17) 


0.077 (28) 


{ad,niv) 


0.10 (91) 


0.073 (94) 


(ci,nid) 


0.058 (12) 


0.057 (13) 


(CI,NIV) 


0.035 (12) 


0.033 (18) 


(nid,niv) 


0.042 (6) 


0.042 (6) 



TABLE VIII: Forecasted uncertainties on isocurvature pa- 
rameters for different cases for the Planck and LSS data 
(Boss and Adept) for fully correlated isocurvature case. The 
percentage improvement in la errors when the LSS data is 
added to the Planck data is shown in brackets. 



tributions and wq. The two degenerate directions are 
almost orthogonal. Here, the inner straight lines repre- 
sent the marginalised error bars obtained by combining 
both Planck and BOSS experiments in the case of a 
cosmological constant. Clearly the ability of LSS data to 
measure isocurvature modes is related to the information 
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FIG. 12: EfTect of combining the CMB and LSS datasets on 
the 1(7 error eUipses for isocurvature contributions and wq. 
We have only represented the isocurvature modes ((ad,niv) 
and (ad,nid)) that determine the main degenerate direction 
with the largest isocurvature fraction. Inner straight lines 
represent the error bars obtained using both CMB and LSS 
experiments, but assuming a cosmological constant. 



provided by the BAO about dark energy. In Figure [131 
we compare the la error ellipses for (wp, (ad,NIv)) and 
(wo, (aDjNID)) that are obtained when we include the 
CMB dataset and add the data from the BOSS redshift 
bins in succession. We see that the BAO data primarily 
serves to reduce the phase space for wq with the largest 
improvement in the wq constraint coming from the sec- 
ond redshift bin. As the redshift increases, the contri- 
bution from dark energy diminishes until matter comes 
to dominate, at which time the impact of dark energy 
on the observables is small. For this reason, the inter- 
mediate redshift bin for Boss centered at z = 0.6 pro- 
vides the best constraints on wq. In figure [Ml we compare 
the sum of the most dominant isocurvature contributions 
^(AD,NIV)^<^°'^^^> +^(AD,NID)P<^°'^^^> to the power 
spectrum at the different redshift bins of the BOSS exper- 
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FIG. 13: Effect of adding different redshift bin datasets on the 
la error ellipses for {wo,Wa), (wo, (ad,niv)) (wq, (ad,nid)). 
Hatched regions on panels (b) and (c) represent the la error 
ellipse for the CMB experiment alone. 



imcnt to their respective error bars. The area between 
the solid (signal) and dotted (error) curves indicates the 
amount of information provided by each bin. Clearly, 
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FIG. 14: Sum of the most dominant isocurvature contribu- 
tions 2<AD,NIV>P<^°'^^^> +^<AD,NID)J'<^°'^'°> to the power 
spectrum for different redsliift bins of the BOSS experiment. 
The blue, green and red solid curves represents respectively 
the isocurvature contribution at redshift z = 0.35, z = 0.6 
and 2 = 3. The dotted lines represent the BOSS error bars for 
different redshift bins. 





AD 








CI ■■■■ 








Bl - - 








NID 








N IV — ■ 























O 2 A 6 S 1 O 



I AD 






^ C 






v'^X Bl - - 






,"\-\ NID -■- 












NV - 



























O 2 4 6 B 1 O 



FIG. 15: Growth function as a function of redshift for all 
modes. Top: on intermediate scales (fc = 0.01 h Mpc~^). 
Bottom: on large scales k = 5.3 x 10~^ h Mpc~^. On small 
and intermediate scales, the growth function is the same for 
all regular modes, while the AD mode grows faster than the 
isocurvature modes on large scales. We have normalized the 
growth to unity at z — 10. 



this combination of isocurvature parameters is best con- 
strained from the measurement of the galaxy power spec- 
trum at z = 0.6 for this particular experiment. Further- 
more, the differing shapes of the signal curves suggests 
that complementary information is available at different 



rcdshifts. Hence, the measurement of the BAO scale at 
different redshifts between decoupling and today helps to 
constrain the isocurvature modes. 

We note that in this study we have assumed exact 
knowledge of the nonlinear shift of the BAO location as 
a function of redshift. 

As an aside, we note that information about the initial 
conditions from LSS data does not stem from differences 
in the growth rates for different modes. Figure [15] shows 
the growth function of the perturbations on intermedi- 
ate (top) and on large scales (bottom). It is clear that 
on very large scales, the isocurvature modes grow more 
slowly than the adiabatic modes. This is expected as 
perturbations which are isocurvature in nature only grow 
when they enter the horizon while adiabatic fluctuations 
grow at all times. However on the scales probed by the 
BAO signal, the isocurvature modes and adiabatic modes 
grow at the same rate. 



IV. CONCLUSIONS 

The first detection of the BAO peak in the galaxy cor- 
relation function measured by SDSS opened the door to 
using the clustering of galaxies on scales of ^ 150 Mpc as 
a cosmic yardstick. By comparing the size of the overden- 
sity of baryons at the epoch of recombination predicted 
from theory and calibrated by the CMB, with its size as it 
appears in the large-scale structure of galaxies today, we 
can study the expansion history of the universe. How- 
ever, in order to succeed in making a precise measure- 
ment of the signal we will need the huge volumes probed 
only by the most recent generation of redshift surveys. 
With such precision we hope to reveal the nature of dark 
energy and probe its time evolution if it exists. 

With claims of constraints on dark energy from BAO 
experiments to the level of a few percent, it becomes 
important to check the assumptions made in the post- 
observational analysis. In this paper, we have revisited 
the assumption of pure adiabatic initial conditions and 
considered the impact of allowing isocurvature-adiabatic 
admixtures on the BAO peak and the implications for 
dark energy studies. We have shown that a combination 
of differences in the baryon growth profile that arises due 
to the presence of isocurvature modes and Silk damp- 
ing change both the shape and position of the BAO 
peak. Non-adiabatic initial conditions leave the sound 
speed unchanged but instead alter the development of 
the acoustic waves in the baryon-photon fiuid prior to 
decoupling which modifies the scale on which the sound 
waves imprint on the baryon distribution. 

The degeneracy between the impact of mixed initial 
conditions and the effect of a dynamical dark energy 
model on the BAO peak weakens the potential con- 
straints on the dark energy parameters forecasted for a 
combined Planck and LSS dataset. We found that the 
admission of more general initial conditions which include 
isocurvature modes and their cross-correlations increases 
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the 95% confidence region in {wo,Wa) space by 50% in 
the case of the Boss experiment and thus the assump- 
tion of adiabaticity can lead to the under-estimation of 
the errors on the dark energy parameters. Furthermore, 
if we assume purely adiabatic initial conditions, we run 
the risk of attributing a shift in the peak away from the 
prediction of '--^ 150 Mpc for a ACDM universe to a non- 
A dark energy model. We have shown that this can lead 
to a bias in the estimates of the dark energy parame- 
ters, leading to a several a incorrect confirmation of A or 
detection of non-A. 

On a positive note, the change in the BAO peak in 
isocurvature models indicates that there is useful infor- 
mation in the galaxy correlation function on the nature of 
the primordial perturbations even when simultaneously 
measuring dark energy equation of state parameters. We 
find that the use of the LSS data in addition to the CMB 
data substantially improves our ability to measure the 
contributions of different modes to the initial conditions. 
The matter power spectrum constrains the dark energy 
parameters and in so doing breaks the degeneracy in the 



isocurvature-dark energy parameter space. Furthermore, 
even when assuming w = —1, the degenerate parame- 
ter combinations in the CMB and LSS are different. A 
similar conclusion is reached in [s^ despite differences in 
their analysis compared with ours, such as their choice 
of non-adiabatic fiducial model, the inclusion of spatial 
curvature and a different LSS experiment. 
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VI. APPENDIX 



The derivatives of the matter power spectrum with re- 
spect to the isocurvature parameters have not been pre- 
sented in the literature before and we show them here for 
the benefit of the reader, in addition to the derivatives 
with respect to the cosmological parameters. 
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FIG. 17: Logarithmic derivatives of P{k) with respect to the isocurvature parameters for different redshifts: z — 0.35 (solid 
black), z = 0.6 (dotted red) and z = 3 (dashed blue). An adiabatic fiducial model is assumed. 
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